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A path integral formulation is developed for the dynamic Casimir effect. It allows us to study 
arbitrary deformations in space and time of the perfectly reflecting (conducting) boundaries of a 
cavity. The mechanical response of the intervening vacuum is calculated to linear order in the 
frequency-wavevector plane. For a single corrugated plate we find a correction to mass at low 
freqüències, and an effective shear viscosity at high freqüències; both anisotropic. For two plates 
there is resonant dissipation for all freqüències greater than the lowest optical mode of the cavity. 
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The Standard Casimir effect Q is a macroscopic man- 
ifestation of quantum fluctuations of vacuum. The mod- 
ified boundary conditions of the electromagnètic field in 
the space between two parallel conducting plates changc 
zero point vacuum fluctuations, resulting in an attractive 
force between the plates, which is within the scnsitivity of 
current experimental force apparatus 0]. Thus, by mea- 
suring the mechanical force between macroscopic bod- 
ies, it is in principle possible to gain information about 
the behavior of the quantum vacuum. Although less well 
known than its static counterpart, the dynamical Casimir 
effect, describing the force and radiation from moving 
mirrors has also garnered much attention This is 

partly due to connections to Hawking and Unruh effects 
(radiation from black holes and accelerating bodies, re- 
spectively), suggesting a deeper link between quantum 
mechanics, relativity, and cosmology |ïü| . 

The creation of photons by moving mirrors was first 
obtained by Moore || for a 1 dimensional cavity. Fulling 
and Davis M demonstrated that there is a corresponding 
force, even for a single mirror, which depends on the third 
time derivative of its displacement. Thcsc computations 
take advantage of conformal symmetries of the 1+1 di- 
mensional space time, and can not be easily gener alized 
to higher dimension. Furthermore, the calculated force 
has causality problems reminiscent of the radiation reac- 
tion forces in classical electron theory @). It has been 
shown that this problem is an artifact of the unphysical 
assumption of perfect reflectivity of the mirror, and can 
be resolved by considering realistic frequency dependent 
reflection and transmission from the mirror ||. 

Another approach to the problem starts with the fluc- 
tuations in the force on a single plate. The fluctuation- 
dissipation theorem is then used to obtain the mechanical 
response function || , whose imaginary part is related to 
the dissipation. This method does not have any causality 
problems, and can also be extended to higher dimensions. 
(The force in 1+3 dimensional space-time depends on the 
fifth powcr of the motional frequency.) The emission of 
photons by a perfect cavity, and the observability of this 



energy, has been studied by different approaches j7p|j). 
The most promising candidate is the resonant production 
of photons when the mirrors vibrate at the optical reso- 
nance frequency of the cavity |ï(Í[j . A review, and more 
extensive references are found in Ref. |ÏÏJ . More recently, 
the radiation due to vacuum fluctuations of a collapsing 
bubble has been proposed as a possible explanation for 
the intriguing phenomenon of sonoluminescense |p^ , ^3| . 

In this letter we present a path integral formulation, 
applicable to all dimensions, for the problem of perfectly 
reflecting mirrors that undergo arbitrary dynamic defor- 
mations |Ï4|] . We calculate the frequency-wavevector de- 
pendent mechanical response function, defined as the ra- 
tio between the induced force and the deformation ficld, 
in the linear regime. From the response function we 
extract a plethora of interesting results, some of which 
we list here for the specific example of lateral vibrations 
of uniaxially corrugated plates: (1) A single plate with 
corrugations of wavenumber k, vibrating at freqüències 
lü <C cfc, obtains anisotropic corrections to its mass. (2) 
For ui 3> ck, there is dissipation due to a frequency de- 
pendent anisotropic shear viscosity. (3) A second plate 
at a separation H modifies the mass renormalization by 
a function of kH, but does not change the dissipation for 
freqüències lo 2 < (ck) 2 + (irc/H) 2 . (4) For all freqüències 
higher than this first optical normal mode of the cavity, 
the mechanical response is infinite, implying that such 
modes can not be excited by any finite external force. 
This is intimately connected to the resonant particle cre- 
ation reported in the literature iQ-fu]]. (5) A phase angle 
9 between two similarly corrugated parallel plates results 
in Josephson-like effects: There is a static force propor- 
tional to sin(#), while a uniform relative velocity results 
in an oscillating force. (6) We calculate a (minute) cor- 
rection to the velocity of capillary waves on the surface 
of mercury due to a small change in its surface tension. 

Our approach is a natural generalization of the path in- 
tegral methods developed by Li and Kardar to study 
fluctuation-induced interactions between deformed man- 
ifolds embedded in a correlated fluid. Such interactions 
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result from thermal fluctuations of the fluid. These meth- 
ods are readily generalized to zero point quantum fluc- 
tuations of a field, taking advantage of the path integral 
quantization formalism. Since in the Euclidian path in- 
tegral formulation the space and time coordinates are 
equivalent, deformations of the boundaries in space and 
time appear on the same footing. As is usual, we simplify 
the problem by considering a scalar field <f> (in place of 
the electromagnètic vector potential |Ï6| ) with the action 

S=^J à d X d^{X)d^{X), (1) 

where summation over fj, = 1, • • • , d is implicit. Following 
a Wick rotation, imaginary time appears as another co- 
ordinate Xd = ict in the ci-dimensional space-time. We 
would like to quantize the field subject to the constraints 
of its vanishing on a set of n manifolds (objects) defined 
by X — X a (y a ), where y a parametrize the ath manifold. 
Following Ref. |ïq ], we implement the constraints using 
delta functions, and write the partition function as 



Z = J Wnn f (^ X «W)) exp|--5Mj. (2) 

The delta functions are next represented by integrals 
over Lagrange multiplier fields. Performing the Gaussian 
integrations over <p(X) then leads to an effective action 
for the Lagrange multipliers which is again Gaussian [ fÏ5| . 
Evaluating Z is thus reduced to calculating the logarithm 
of the determinant of a kernel. Since the Lagrange mul- 
tipliers are defined on a set of manifolds with nontrivial 
geometry, this calculation is generally complicated. To 
be specific, we focus on two parallel 2d plates embedded 
in 3+1 space-time, and separated by an average distance 
H along the a^-direction. Deformations of the plates 
are parametrized by the height functions /ii(x, í) and 
/i2(x, í), where x = [x\, x-z) denotes the two lateral space 
coordinates while t is the time variable. Following Ref. 
]Ï5| , ln Z is calculated by a perturbative series in powers 
of the height functions. The resulting expression for the 
effective action (in real time), defined by S e s = —ih\nZ, 
after eliminating h independent terms, is 

(3) 



-A-(q,w) (h 1 (q,üj)h 2 (-q, -üj) + h^-q, -uj)h 2 (q,uj))} + 0(h 3 ). 



The kernels A±(q,u), that are closely related to the 
mechanical response of the system (see below) , are func- 
tions of the separation H, but depend on q and to only 
through the combination Q 2 = q 2 — uj 2 /c 2 . The closed 
forms for these kernels involve cumbersome integrals, and 
are not very illuminating. Instead of exhibiting these 
formulas, we shall describe their behavior in various re- 
gions of the parameter space. In the limit H — » oo, 
A°°(q,iü) = 0, and 

í ~ 360^? _ w2 ) 5/2 fOT U < C «' 

A°?{q,u) = { (4) 

l l K^(^ 2 - c V) 5/2 forocg, 

where sgn(w) is the sign function. While the effective ac- 
tion is real for Q 2 > 0, it becomes purely imaginary for 
Q 2 < 0. The latter signifies dissipation of energy ||, pre- 
sumably by generation of photons || . It agrees precisely 
with the results obtained previously H for the special 
case of flat mirrors (q = 0). (Note that dissipation is 
already present for a singlo mirror.) 

In the presence of a second plate (i.e. for finite H), the 
parameter space of the kernels subdivides into three dif- 
ferent regions as depicted in Fig. 1. In region I (Q 2 > 
for any H), the kernels are finite and real, and hence there 



is no dissipation. In region lla where —ir 2 /H 2 < Q 2 < 0, 
the í/-independent part of A + is imaginary while the 
iï-dependent parts of both kernels are real and finite. 
(This is also the case at the boundary Q 2 = —it 2 /H 2 .) 
The dissipation in this regime is simply the sum of what 
would have been observed if the individual plates were 
decoupled, and unrelated to the separation H . By con- 
trast, in region Ilb where Q 2 < —ir 2 /H 2 , both kernels 
diverge with infinite real and imaginary parts (ÏtJ . This 
iï-dependent divergence extends all the way to the nega- 
tive Q 2 axis, where it is switched off by a Í/H 5 prefactor. 

As a concrete example, let us examine the lateral 
vibration of plates with fixed roughness, such as two 
corrugated mirrors. The motion of the plates enters 
through the time dependences /ii(x, í) = hi(x — r(í)) 
and /i2(x, t) = /i2(x); i.e. the first plate undergoes lat- 
eral motion described by r(í), while the second plate is 
stationary. The lateral force exerted on the first plate 
is obtained from /i(í) = SS e g/5ri(t). Within linear re- 
sponse, it is given by 

/ i (w)=x«Hr i (a;) + / i («) ) (5) 

where the "mechanical response tensor" is 



Xij(u) = te I qiqA[A + (q,u;) - A + (q,0)] ^(q)] 2 + ±A.(q, 0) (^i(q)^(-q) + M-q)te(q))} - (6) 

and there is a residual force 

fir í d 2 n 

/°H = -— 2nS(u J ) J — ^ i ft A_(ç,0) (/n(q)fta(-q) - ^(-^^(q)) . (7) 
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For a single corrugated plate with a deformation 
h{x) — (fcosk-x, we can easily calculate the re- 
sponse tensor using the explícit formulas in Eq.(||). In 
the limit of w <C ck, expanding the result in pow- 
ers of lü gives Xij — bm^uí 2 + 0(cu 4 ), where Sniij — 
ATid 2 k 3 L·iL·j /(2887r 2 c), can be regarded as corrections to 
the mass of the plate. (Cut-off dependent mass correc- 
tions also appear, as in Ref. |ÏI}| .) Note that these mass 
corrections are anisotropic with 8m\\ = Afrk 5 d 2 / {288ir 2 c) 
and Sm± — 0. Parallel and perpendicular components 
are defined with respect to k, and A denotes the area of 
the plates. The mass correction is inherently very small: 
For a macroscopic sample with d ~ X = 2n/k w lmm, 
density ~ 15gr/cm , and thickness t w 1 mm, we find 
Sm/m ~ 10 -34 . Even for deformations of a micro- 
scopic sample of atòmic dimensions (close to the limits 
of the applicability of our contimmm representations of 
the boundaries), Sm/m can only be reduced to around 
10~ 10 . While the actual changes in mass are immeasur- 
ably small, the hope is that its anisotropy may be more 
accessible, say by comparing oscillation freqüències of a 
plate in two orthogonal directions. 

For lü 3> ck the response function is imaginary, and we 
define a frequeney dependent cffective shear viscosity by 
Xijiio) = —iujrjij(iü). This viscosity is also anisotropic, 
with 77|| {lü) = hAk 2 d 2 Lü 4 /{720Tr 2 c 4 ), and tí_l(w) = 0. 
Note that the dissipation is proportional to the fifth time 
derivative of displacement, and there is no dissipation 
for a uniformly accelerating plate. However, a freely 
oscillating plate will undergo a damping of its motion. 
The characteristic decay time for a plate of mass M is 
r k 2M/r\. For the macroscopic plate of the previous 
paragraph, vibrating at a frequeney of lü 2ck (in the 
10 12 Hz range), the decay time is enormous, r ~ 10 18 s. 
However, since the decay time seales as the fifth power 
of the dimension, it can be reduced to 10~ 12 s, for plates 
of order of 10 àtoms. However, the required freqüències 
in this case (in the 10 18 Hz range) are very large. Also 
note that for the linearized forms to remain vàlid in this 
high frequeney regime, we must require very small ampli- 
tudes, so that the typical velocities involved v ~ r^oj, are 
smaller than the speed of light. These difhculties can be 
somewhat overcome by considering resonant dissipation 
in the presence of a second plate. 

With two plates at an average distance H, the results 
are qualitatively the same for freqüències less than the 
natural resonance of the resulting cavity. There is a 
renormalization of mass in region I, and dissipation ap- 
pears in region lla, of Fig. 1. However, the mass renor- 
malization at low freqüències {lü <C ck) is now a func- 
tion of both k and H, with a crossover from the single 
plate behavior for kH ~ 1. In the limit of kH <C 1, 
we obtain 5m\\ = hABk 2 d 2 / 48cH 3 and Sm± = 0, with 
B = —0.453. Compared to the single plate, there is an 
enhancement by a factor of {kH)~ 3 in Sm\\ . The effective 
dissipation in region lla is simply the sum of those due 
to individual plates, and contains no H dependence. 



There are additional interesting phenomena resulting 
from resonances. We find that both real and imaginary 
parts of A±{q,Lü), diverge for lü 2 /c 2 > q 2 + n 2 /H 2 . In 
the example of corrugated plates, we replace q by k to 
obtain a continuous spectrum of freqüències with diverg- 
ing dissipation. Related effects have been reported in 
the literature for 1+1 dimensions ]7| [ÏOj] , but oceuring 
at a discrete set of freqüències uj n = nirc/H with integer 
n > 2. These resonances oceur when the frequeney of the 
external perturbation matches the natural normal modes 
of the cavity, thus exciting quanta of such modes. In one 
space dimension, such modes are characterized by a dis- 
crete set of wavevectors that are integer múltiples oiir/H. 
The restriction to n > 2 is a consequence of quantum 
electrodynamics being a 'free' theory (quadratic action): 
only two-photon states can be excited subject to conser- 
vation of energy. Thus the sum of the freqüències of the 
two photons should add up to the external frequeney |j) . 

In higher dimensions, the appropriate parameter is the 
combination lü 2 j c 2 — q 2 . From the perspective of the ex- 
cited photons, conservation of momentum requires that 
their two momenta add up to q, while energy conser- 
vation restriets the sum of their freqüències to lü. The 
in-plane momentum q, introduce a continuous degree of 
freedom: the resonance condition can now be satisfied 
for a continuous spectrum, in analogy with optical res- 
onators. In Ref. [9], the lowest resonance frequeney is 
found to be 2ixc/B which seems to contradict our pre- 
diction. However, the absence of lü\ = ttc/H in 1+1 D 
is due to a vanishing prefactor which is also present 
in our calculations. However, in exploring the continu- 
ous frequeney spectrum in higher dimensions, this single 
point is easily bypassed, and there is a divergence for 
all freqüències satisfying lo 2 jc 2 > q 2 +ir 2 /H 2 , where the 
inequality holds in its strict sense. 

Resonant dissipation has profound consequences for 
motion of plates. It implies that due to quantum fluctu- 
ations of vacuum, components of motion with freqüències 
in the range of divergences cannot be generated by any 
finite external forcel The imaginary parts of the kernels 
are proportional to the total number of excited photons 
||. Exciting these degrees of motion must be accompa- 
nied by the generation of an infinite number of photons; 
requiring an infinite amount of energy, and thus impossi- 
ble. However, as pointed out in Ref. ||], the divergence is 
rounded off by assuming finite reflectivity and transmis- 
sivity for the mirrors. Hence, in practice, the restriction 
is softened and controlled by the degree of ideality of the 
mirrors in the frequeney region of interest. 

We shall next examine the constant term in Eq.(Q). 
For two plates corrugated at the same wavelength, with 
deformations /ii(x) = d\ cos(k • x) and /i2(x) = c?2 cos(k • 
x + a), there is a (time independent) lateral force 



hcA 

F dc = ( fc > 0)krfirf 2 sin a, 



(8) 



which tends to keep the plates 180 degrees out of phase, 
i.e. mirror symmetric with respect to their mid-planc. 
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The dependence on the sine of the phase mismatch is 
reminiscent of to the DC Josephson current in supercon- 
ductor junctions, the force playing a role analogous to 
the current in SIS junctions. There is also an analog for 
the AC Josephson effcct, with vclocity (the variable con- 
jugate to force) playing the role of voltage: Consider two 
corrugated plates separated at a distance H, described by 
hi(x,t) = di cos[k-(x — r(í))] and /i2(x, í) = di cos[k-x]. 
The resulting force at a constant velocity (r(í) = v t), 

hcA 

F ac = — A_(M)kdid 2 sin[(k·v) *], (9) 

oscillates at a frequcncy u) = k • v. Actually both effects 
are a consequence of the attractive nature of the Casimir 
force. It would be difficult to separate them from similar 
forces resulting from say, van der Waals attractions. 

As a final example, we study the capillary waves on 
the surface of mercury, with a conducting plate placed at 
a separation H above the surface . The low frequency- 
wavevector expansion of the kernel due to quantum fluc- 
tuations in the intervening vacuum, starts with quadratic 
forms q 2 and u> 2 . These terms result in corrections to 
the (surface) mass density by 5p = hB/4S,cH 3 , and to 
the surface tension by Ser = hcB/48H 3 . The latter cor- 
rection is larger by a factor of (c/c s ) 2 , and changes the 
velocity c s , of capillary waves by <5c s /c° = HcB/96aH 3 , 
where a is the bare surface tension of mercury. Taking 
H ~ lmm and a <~ 500 dynes/cm, we find another very 
samll correction of Sc s /c® ~ 10~ 19 . 

In summary, we have developed a path integral for- 
mulation for the study of quantum fluctuations in a cav- 
ity with dynamically deforming boundaries. As opposed 
to previous emphasis on spectra of emitted radiation, 
we focus on the mechanical response of the vacuum. 
Most of the predicted phenomena, whilc quite intrigu- 
ing theoretically, appear to be beyond the reach of cur- 
rent experiment: the most promising candidates are the 
anisotropy in mass, and resonant dissipation. The path 
integral method is quite versatile, and future extensions 
could focus on non-linear response, other geometries (e.g. 
wires), the gauged electromagnètic field, and calculations 
of emitted spectra using correlation functions. 
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